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ABSTRACT 


The  Wavenumber  Technique  (WT)  is  a  relatively  new  method 
of  underwater  sound  transmission  analysis.  One  aspect, 
source  depth  determination,  is  studied  to  evaluate  its 
validity  and  test  environmental  and  acoustic  sensitivity. 
The  horizontal  wavenumber  spectrum  is  analyzed  to  determine 
null  spacings  in  wavenumber  space,  which  indicates  source 
depth  by  the  Lloyd's  Mirror  interference  effect.  Comparison 
of  this  theory  with  cases  of  an  isospeed  sound  profile, 
fully  absorbing  bottom,  and  flat  totally-reflecting  surface 
shows  excellent  agreement  for  several  parametric  variations. 
Cases  with  realistic  sound  speed  profiles  and  partially 
absorbing  bottoms  generally  agree  with  theory,  but  a  dis¬ 
tinct  bias  is  observed.  Source  depth  determination  curves, 
which  relate  the  scaled  wavenumber  spectral  intensity  null 
spacing  to  the  source  depth,  are  presented  for  comparison 
with  theory  and  sensitivity  analysis.  in  example  is  given 
for  suggested  application  of  source  depth  determination. 
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1.  IHTRODOCTION 


The  Havenuaber  Technique  (HT>  is  a  new  method  for  the 
solution  of  several  problems  related  to  the  analysis  of 
underwater  sound  transmission  (Lauer,  1979)  .  Its  basis  is  the 
analysis  of  acoustic  propagation  in  the  wavenumber  domain, 
which  has  been  described  by  OiNapoli  (1971)  and  was  an 
intermediate  step  in  the  computation  of  transmission  loss 
versus  range  in  the  Fast  Field  Prograi  (FFP) .  DiNapoli 
(1977)  later  used  the  wavenumber  domain  to  study  the  impe¬ 
dance  of  the  ocean-bottom  interface.  Lauer  (1979)  first 
introduced  the  WT  and  its  applications  to  passive  localiza¬ 
tion  and  tracking  and  multipath  decomposition.  He  developed 
a  range  prediction  curve  based  upon  peak  wavenumbers  in  the 
transform  spectrum  fcr  successive  range  increments  from  an 
entire  pressure  range  field. 

These  investigations  presented  highly  beneficial  and 
promising  results  from  the  wavenumber  depicrion,  which  indi¬ 
cates  that  additional  research  is  appropriate.  One  particu¬ 
lar  aspect  of  the  WT  -  determination  of  source  depth  -  will 
be  investigaxed  here  in  order  to  evaluate  its  validity  in 
comparison  with  the  Lloyd's  Hirror  effect.  Its  relation¬ 
ships  to  oceanic  and  acoustic  variations  will  also  be 
examined . 
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The  waeenaebec  technique  as  first  described  by  Lauer 
(1979)  vas  based  upon  input  froa  the  FFP,  but  it  has  been 
adapted  for  use  with  the  Parabolic  Equation  (PE)  method  of 
propagation  loss  deteraination.  The  PB-based  WT  is  currently 
being  investigated  by  Lauer  and  others  at  the  Naval  Ocean 
Research  and  Develops ent  Activity  (NORDA)  and  was  the  aodel 
used  in  this  research  at  tha  Naval  Postgraduate  School 
(NPS) .  Descriptions  of  the  PB  and  HT  will  be  presented  to 
display  their  sensitivities  to  variations  of  several  aodel 
parameters.  Also,  expected  relationships  of  the  WT  to  other 
current  propagation  prediction  models  will  be  discussed  to 
consider  how  results  might  vary  based  upon  different  inputs. 
Particular  attention  will  be  given  to  algorithm  differences, 
degrees  of  approximation,  and  treatment  of  the  ocean  bottom. 

The  PE-based  WT  will  then  be  exercised  for  several 
parameter  variations  to  evaluate  qualitatively  the  consis¬ 
tency  of  the  technique.  This  study  will  also  demonstrate 
the  response  of  the  PE-based  WT  to  geometric  and  actual 
oceanic  conditicns.  Prediction  curves  for  source  depth 
determination  based  upon  these  variations  will  be  presented 
and  discussed. 
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II.  iA7EN0HBgR  TECHNIQPE 


A.  GENERAL  DESCRIPTION 

The  NT  is  a  natoral  alternative  for  signal  depiction 
because  it  provides  inforoation  on  the  directions  of  energy 
arrivals.  It  allows  reasonably  direct  physical  interpreta¬ 
tion  of  ocean  acoustic  processes  and  their  relationships  to 
environaental  conditions  (Lauer,  1979)  .  Lauer  further  states 
that  the  wavenunber  spectrue  plot  is  clearer  and  more 
inforaative  than  typical  curves  of  propagation  loss  versus 
range.  The  wavenumber  K  relates  the  angular  frequency  and 
*he  sound  speed  by 


K  =  w  =  2Trf 
C  C 


and  it  has  horizontal  and  vertical  components 


K  =  k  +  k  , 
r  2 


The  reference  wavenuaber  is  determined  by  the  sound  speed 


mini  nun 


K  =0) 
0  — 
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The  WT  will  first  be  described  without  reference  to  any 
particular  propagation  model  to  present  the  general  nature 
of  the  algorithm.  Fig.  1  is  the  flow  iiagram  for  the  WT. 
The  left  and  right  sides  of  this  flow  diagram  correspond  to 
the  intensity  and  wavenumber  axes,  respectively,  of  the 
spectrum  plot  that  will  result  from  the  HT.  This  figure  and 
the  following  discussion  were  developed  from  Lauer  (1979  and 
1982a)  and  augmented  by  several  telephone  conversations  with 
B.  Lauer,  T.  Lawrence,  and  R.  Evans  of  NORDA  between  July 
and  December,  1982. 

The  WT  calculates  spectral  intensity  from  the  real  and 
imaginary  parts  of  the  acoustic  pressure 

p  =  ^(R,Z)  . 


The  complex  psi  versus  range  field  at  a  specified  depth  is 
the  propagation  model  product  which  determines  the  pressure. 
The  complex  pressure  field  can  be  modified  for  effects  such 
as  attenuation,  depending  upon  whether  idealized  or  realise 
tic  conditions  are  being  investigated.  A  Fast  Fourier 
Transform  (FFT)  is  then  applied  to  yield  transformed  complex 
pressure.  The  spectral  intensity  is  de-^ermined  at  each 


2D 


r 


Figure 


WAVENUMER  TECHNIQUE  FLOW  D2ACRAM 


.  Wavenumber  Technique  Flow  Oiagrio. 
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wavenuaber  increaent  from  the  transformed  complex  pressure 
field  by 

Is  *  "  I^'t'  • 


This  intensity  field  is  scaled  arbitrarily  for  satisfactory 
display  of  the  relative  maxima  and  minima  to  define  the 
intensity  nulls. 

Additional  output  parameters  from  the  PE  which  are 
required  for  the  HT  include  the  range  step,  the  number  of 
range  points,  frequency,  minimum  sound  speed  and  field 
depth.  These  are  the  selected  variables  noted  in  figure  1. 
The  range  step  is  selected  to  provide  sufficient  resolution 
to  determine  the  wavenumber  null  spacing  on  the  spectrum 
plots.  The  initial  radial  component  of  the  wavenumber  for 
the  spectrum  is  calculated  from  DiNapoli  (1971)  by 

k  :  K  -  2Tr 
ro  o  — 

AR 


where 


K 

0 


ro 


+  k. 
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A 

is  will  be  seen  in  section  II-B,  if  a  "scaled”  wavenua- 
ber  beta  is  used  for  the  horizontal  axis,  analysis  is  sia- 
plified  since  the  nulls  in  the  spectra  should  be  more  evenly 
spaced  (Lauer,  1979) .  The  physical  meaning  of  beta  will  be 
described  in  the  next  section.  The  resulting  scaled  wave¬ 
number  spectra  are  then  plotted.  The  wavenumber  null  spacing 
for  each  increaental  source  depth  is  determined  from  each 
scaled  spectrum.  The  source  depth  determination  curve  is 
generated  by  plotting  these  null  spacings  as  a  function  of 
source  depth,  as  suggested  by  Lauer  (1979)  .  The  null  spac¬ 
ing  of  an  analyzed  received  signal  could  be  compared  with 
this  determination  curve  to  infer  the  source  depth  of  the 
signal.  This  application  will  be  discussed  later.  The 
source  depth  determination  curve  will  be  specific  for  the 
acoustic  and  oceanic  descriptions  of  the  medium. 

B.  PHYSICAL  DESCRIPTION 

The  preceding  description  was  intended  to  provide  the 
overall  concept  of  the  NT  whereas  whar  follows  will  explain 
the  physical  basis  for  the  technique.  The  underlying  princi¬ 
ple  of  the  WT  as  described  by  Lauer  (1982a)  is  the  Lloyd's 
Mirror  affect,  which  gives  the  acoustic  field  by 

p  =  E(R)  ^ 
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where 


E(R)  =  exp  {iK  /  -  Zg)^}  -  exp  { iK  /  R^  +  )  ^ } 

/  r2  +  (Zj^  -  Zg)^  /  R2  +  (Zj^  + 

for  a  point  source  in  a  semi- infinite  ledium  of  constant 
sound  speed  with  a  pressure  ralease  surface.  This  is  the 
Lloyd's  Mirror  field  in  range  space,  whereas  the  Lloyd's 
Mirror  field  in  wavenumber  space,  applicable  to  this 
research,  is  given  by  Lauer  (19  82)  as 

F(K)  =  sinCBZg)  exp(iBZj^)  , 

8 


where  beta  is  defined  by 

3  =  (K2_k2)''S 

or* 


and  E(R)  and  F(K)  are  related  by  the  3essel  transform  pair 


E(R)  =  /"2  FCO  J  (KR)  K  dK 
0  o 


and 

F(K)  =  1  /"“eCR)  J  (KR)  R  dR 
o  o 
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For  this  research  the  Bessel  function  vas  approxioat ad  as  a 
trigonoaetric  function  so  that  the  FFT  could  be  utilized. 
This  is  considered  acceptable  since  tha  source  and  receiver 
are  greatly  separated,  according  to  Lauer  (1982b).  Inspec¬ 
tion  of  f(K)  reveals  that  the  nulls  of  the  spectrum  will  be 
equally  spaced  at  intervals  of 

AS  =  IV  ^ 

"s 

This  equal  spacing  in  beta  is  the  )cey  element  of  the  appli¬ 
cation  of  the  WT  to  source  depth  determination  (Lauer , 1 979) . 

Only  the  direct  and  surface -ref lectai  waves  will  inter¬ 
act  at  a  receiver  location  in  a  simple  illustration  of  the 
Lloyd's  Wiinror  effect  with  a  smooth  surface.  This  interac¬ 
tion  will  be  destructive  or  constructive,  as  a  result  of 
their  relative  phases  (Kinsler  at  al  ,1982)  and  will  produce 
nulls  or  peaks,  respectively,  in  the  wavenumber  spectrum. 
Fig.  2  depicts  the  Llcyd's  Mirror  effect  geometry  which  will 
be  used  to  study  tesr  case  outputs  for  null  spacing  and  sub¬ 
sequent  source  depth  determination.  The  surface-reflected 
(SH)  wave  travel  distance  will  increase  at  a  "faster  rate 
than  tha  direct  (D)  wave  distance  as  the  source  is  moved 
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Surface 


BOttOM 


Surface 


Bottoa 


deeper  in  the  water  cclunn  (d)  while  'rhe  receiver  depth  (h) 
is  stationary.  This  leads  to  different  arrival  patterns  and 
nulls  at  acre  wavenuabers.  Thus  the  general  pattern  of 
decreased  spacing  in  beta  with  increasing  source  depth  is 
observed . 

Bottom  interacting  waves  will  be  neglected  for  idealized 
cases.  This  will  be  accomplished  by  tie  use  of  a  fully 
absorbing  bottom.  These  waves  will  be  important,  however, 
in  shallow  water  or  at  greater  ranges,  k  few  test  cases  will 
have  realistic  bottoms.  Various  propagation  models  handle 
the  bottom  in  different  ways,  so  the  inclusion  of  a  bottom 
will  be  discussed  later.  Cases  in  shallow  water  (less  than 
5000  feet)  or  at  extended  ranges  (greater  than  100  miles) 
will  not  be  considered  here. 

Lauer  (1979)  observed  that  the  wavenumber  specrrum  shews 
a  series  of  modulating  envelopes  which  enclose  spikes  rhat 
are  the  eigenvalues  associated  with  the  normal  modes.  Com¬ 
parison  of  the  wavenumber  spectrum  with  typical  transmission 
loss  curves  reveals  that  the  wavenumber  depiction  is  mere 
coherently  structured.  This  leads  to  a  better  physical 
understanding  of  signal  transmission  because  of  the  rela¬ 
tionship  between  horizontal  wavenumber  and  ray  nheory.  He 
gives  this  as 

^  cos  0  =  2  Trf  COS0  . 
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The  horizontal  wavenuDber  decreases  as  theta  increases.  Thus 


r 


the  wavenuaber  axis  of  the  spectrum  plot  will  include  vave- 
nuabers  associated  with  saaller  elevation  (or  depression) 
angles  in  increasing  order  of  wavenumber  (Lauer,1979) .  It 
is  then  expected,  for  the  case  of  a  neglected  bottom,  thar 
spectrum  amplitude  peaks  could  be  observed  that  would  corre¬ 
spond  to  the  direct  and  surface  reflected  waves. 

C.  iJAVENOMBEE  TECHNIQUE  APPLICATION 

Lauer  (1979)  recommends  a  receiver  that  is  a  single 
omnidirectional  hydrophone  and  a  source-generated  continuous 
wave  (CH)  received  signal  for  source  depth  determination. 
He  further  states  that  in  actual  use  the  single  hydrophone 
could  be  replaced  by  an  array  to  improve  the  sianal-co-noise 
ratio  and  to  obtain  bearing  and  bearing  rate  information. 

An  example  of  a  possible  application  of  the  WT,  source 
depth  determination,  is  presented  in  Fig.  3.  A  propagation 
model  could  be  run  for  a  series  of  incremental  source  depths 
to  yield  the  null  spacing  relationship,  as  was  shown  in  Fig. 
1.  That  curve  would  be  specific  to  the  properties  of  the 
ocean  and  values  for  parameters  such  as  frequency  and 
receiver  depth.  A  received  signal  as  shown  on  the  righr  side 
of  Fig.  3  could  then  be  analyzed.  Its  null  spacing  could  be 
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coaparei  to  the  model-produced  curve  to  infer  the  depth  of 
the  source. 

D.  TECHNIQUE  LIMITATIONS 

The  success  of  the  NT  is  related  directly  to  the  accu¬ 
racy  of  the  predicted  pressure  field  produced  by  the  propa¬ 
gation  model.  There  are  several  possible  differences 
between  propagation  models  which  might  affect  the  pressure 
field.  Some  of  the  more  significant  include  approximations 
of  the  wave  equation,  bottom  affects,  initial  source  func¬ 
tions  and  input  limitations.  These  will  be  compared  for  sev¬ 
eral  models  to  determine  relevant  restrictions.  The  pressure 
field  will  also  respond  to  changes  in  the  ocean.  A  series 
of  different  scenarios  will  be  run  to  provide  an  initial 
estimate  of  the  sensitivities  to  some  of  these  factors. 

Any  additional  limitations  will  be  related  to  computer 
processing  time  and  storage  requirements.  These  factors  will 
not  be  addressed  directly  in  this  research  because  the  pro¬ 
cessing  following  the  propagation  model  is  minimal.  There  is 
also  no  requirement  tc  store  the  complex  pressure  fields  or 
their  transforms  once  the  null  spacing  has  been  determined. 
Thus  the  NT  computer  utilization  is  directly  proportional  to 
the  time  and  storage  required  for  the  propagation  models. 
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III.  PBOPAGATION  MODELS 

A.  INTSODDCTION 

Modeling  the  propagation  of  sound  in  the  ocean  is  con- 
plicated  by  the  variability  of  the  ocean,  the  greet  range  of 
freguencies  of  interest,  and  the  many  applications  of  sound 
propagation  (DiUapoli  and  Deaven port ,  197 9|  .  it  is  not  sur¬ 
prising  that  no  single  nodal  is  presently  capable  of 
addressing  all  of  these  variations;  there  are  many  semi-re¬ 
strictive  methods.  In  general  all  models  consider  the  sound 
speed  to  be  a  function  of  the  spatial  coordinates  and  inde¬ 
pendent  of  time  (DiNapoli  and  0 ea venport, 1979) .  Only  the 
frequency,  water  column  geometry  and  ooeanic  description 
will  affect  propagation. 

All  models,  regardless  of  application,  can  be  divided 
into  two  classes:  range  independent  or  range  dependent. 
Bangs  independent  models  assume  that  the  ocean  is  cylindri- 
cally  symmetrical,  the  speed  of  sound  is  a  function  only  of 
depth  and  all  boundaries  are  planar  perpendicular  to  the 
depth  axis  (DiNapoli  and  Deavenport,  1979) .  Bange  dependent 
models  can  provide  better  approximations  of  real  condixions 
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by  allowing  the  sound  speed  to  be  a  function  of  two  or  three 
spatial  variables.  They  also  do  not  require  planar  bounda¬ 
ries  perpendicular  tc  the  depth  axis,  thus  cylindrically 
syaaetric  bottoa  topography  can  be  included  (DiMapoli  and 
Deavenport , 1979) .  A  range-dependent  model,  the  split-step 
Fast  Fourier  Transform  PE  (SSPFT)  was  the  numerical  model 
used  in  this  research  to  generate  the  pressure  fields  for 
the  wavenumber  analyses.  This  PE  will  be  qualitatively  com¬ 
pared  to  two  other  range-dependent  algorithms,  a  range-modi¬ 
fied  Normal  Mode  (NM)  model  and  a  Finite  Difference  (FD) 
model.  It  will  alsc  be  compared  tc  a  range-independent 
method,  the  Fast  Field  Program  (FFP)  model.  The  range-de- 
pendent  capability  cf  the  SSFFT  was  not  used  in  this 
research  since  this  was  an  initial  investigation.  Its  com¬ 
parison  with  both  range  dependent  and  independent  models 
may,  however,  indicate  possible  model  effects  on  the  MT 
which  could  be  studied  later. 

In  general  all  models  attempt  to  solve  the  acoustic  wave 
equation.  This  equation  is  valid  for  sound  propagation  in 
fluid  media  with  specified  boundary  conditions 
(0riclc,1 975,p.  1 14)  .  It  becomes  the  reduced  elliptic  Helm¬ 
holtz  equation  in  spatial  coordinates  under  the  assumption 
of  harmonic  time  dependence. 
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B.  P&B&BOUC  EQOiTION 
1 •  Introduction 

The  PE  algoritha  used  la  this  research  was  developed 
by  Brock  (1978)  .  The  basic  feature  of  this  PE  is  the 
replaceaent  of  the  re  diced  elliptic  Relaholtz  equation  by  a 
parabolic  partial  differential  equation  that  utilizes  the 
Tappert- Hardin  split- step  Fourier  algoritha  (SSPFT)  for  the 
nuaerical  integration  ( Brocks  19 78) .  The  original  computer 
code  as  described  by  Brock  has  been  aodified  to  incorporate 
the  source  and  attenuation  descriptions  ot  Tatro  (1977)  and 
to  include  an  ocean  bottoa  by  Stieglitz  ei  al  (1979).  This 
is  the  version  of  the  PE  that  is  currently  available  at  NPS 
and  was  the  code  (PEHODEL)  used  for  this  research.  (It  must 
be  emphasized  that  some  difficulty  was  encountered  in  trying 
to  ascertain  the  correct  description  of  the  NPS  version 
because  of  inadequate  and  erroneous  comment  card  documenta¬ 
tion  in  the  source  code  and  incomplete  manuscript  descrip¬ 
tions.) 

This  sec*:ion  will  discuss  the  SSFFT,  the  original 
Brock  algorithm,  modifications  to  that  code,  and  subsequent 
evaluations  of  the  SSFFT.  Particular  at-.ention  will  be 
given  to  those  aspects  which  ace  expected  to  affect  the  WT. 


2.  Parabolic  Approxiaatlon 


The  parabolic  approximation  to  tha  wave  equation  was 
first  introduced  by  Leontovich  and  Fock  (19U6)  as  a  solution 
for  the  propagation  of  electromagnetic  waves  along  the  sur¬ 
face  of  the  earth.  Its  first  introduction  to  underwater 
acoustics  was  by  Hardin  and  Tappert  (1973)  who  used  the 
split-step  technique  in  conjunction  with  the  FFT  to  solve 
this  particular  apprcximation  to  the  wave  equation.  The 
result  was  the  SSFFT.  Tappert  (1977)  related  the  use  of  the 
parabolic  approximaticn  to  sound  channel  propagation  in  a 
waveguide  that  is  thin  vertically  and  elongated  horizontally 
to  the  range  of  the  first  convergence  zone  or  farther. 

Long-range  propagation,  with  which  the  SSFFT  is 
basically  concerned,  will  usually  be  for  low  frequencies 
because  volume  absorption  increases  strongly  with  frequency 
(Tappert , 1 977) .  The  maximum  elevation  (and  depression) 
angles  of  propagation  for  long  range  transmission  must  be 
small  to  satisfy  the  EE  (Tapper t,  1977)  .  The  derivation  of 
the  SSFFT  will  be  discussed  later  in  a  simplified  version  as 
presented  by  Brock  (1  978). 

It  appears  that  the  SSFFT  is  intended  for  use  in 
situations  where  bottom  interaction  and  surface  scattering 
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would  be  negligible.  Thus  the  use  of  the  SSFFT  to  examine 
the  Lloyd's  Mirror  effect  in  the  real  world  may  be  inappro¬ 
priate  if  the  ocean  surface  is  sufficiently  rough;  addi¬ 
tional  errors  may  occur  in  situations  where  bottom  effects 
are  significant. 

3.  Stasis  Algorithm 
a.  Description 

The  SSFFT  is  valid  for  acoustic  pressure  in  a 
medium  of  constant  density  with  a  monof requency  source  and 
cylindrical  symmetry  about  the  depth  axis  (Brock,  1978)  .  The 
acoustic  pressure  can  be  written  as 

p(R,Z)  =  >p(R,Z)  H  R)  , 

where  the  zero-order  Hankel  function  of  the  first  kind 
relates  the  acoustic  pressure  to  an  outward  propagating 
cylindrical  wave  envelope  function  (Brock,  1 978) .  This  is 
valid  because  of  the  assumption  that  at  low  frequencies  all 
significant  energy  will  propagate  approximately  horizontally 
away  from  the  source  (Tappert,19  77) .  The  asymptotic  form  of 
the  Bankel  function  is 

H  ^(K  R)  »  {2/-nK  R}^  exotiCK  R-ir/4)}  ,KR>>1, 

0  0  0  -  0  0  * 
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if  the  receiver  is  nany  waveleogths  from  the  source 
leads  to  (Brock, 1978) 


This 
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which  is  the  fora  used  by  Brocic  (1978). 
b.  lapleeent  ation 

Brock  (1978)  used  the  SSFFT  to  solve  the  PE 
because  it  has  several  significant  advantages  which  outweigh 
its  disadvantages.  Advantages  include  exponential  accuracy 
in  depth,  second-order  accuracy  in  range,  energy  conserva¬ 
tion,  unconditional  stability  and  computational  efficiency. 
Disadvantages  are  a  uniform  mesh  and  periodic  boundary  con¬ 
ditions  to  satisfy  the  FFT,  and  filtering  of  the  sound  speed 
profile  by  smoothing  discontinuities  to  avoid  spurious  high 
angular- frequency  components  (Br ock,  1 973)  .  Additionally  the 
algorithm  assumes  a  flat  pressure-release  surface,  a  vanish¬ 
ing  field  at  the  maximum  depth  and  a  pseudo-radiation  condi¬ 
tion  at  the  water-bottom  interface  by  saoothly  attenuating 
the  field  (Brock ,  1  978 ) . 

A  numerical  algorithm  of 

=  i(A+B)\p  , 

ip(R+AR,Z)  =  gi-iRO+A)  ^ 


leads  to 

■i;LR+AR,2)  =  FFT~^  {  FFT  (ij^CR,!))  } 


or 


iARK  (n2-l)/2  iARK2/2K 

*(R+AR,2)  =  e  °  FFT"l{e  °FFT(  ij;(R,Z) )  } 


to  facilitate  the  computations  since  a  spatial  FFT  and  its 
inverse  are  required  to  transform  between  depth  and  wavenum¬ 
ber  (Brock,  1 978)  ,  This  procedure  is  implemented  by  two 
alternating  steps,  the  first  of  which  considers  propagation 
in  a  homogeneous  medium  to  account  for  diffraction  and  the 
second  to  account  for  refraction  (Brock, 1 978) .  in  his 
development  Brock  gives  an  alternate  calculation  of  the 
envelope  function  stemming  from 


■^(r+ar,z) 


^iARA/2  ^iARB  iARA/2 


ipC  R ,  Z )  < 


This  differs  slightly  from  the  previous  expression  contain¬ 
ing  the  FFT  and  inverse  FFT.  Thus  it  is  not  currently  known 
which  relationship  is  used  by  Brock;  this  could  be  resolved 
later  by  a  detailed  examination  of  the  compu-er  coda.  In 
any  event  as  McDaniel  (1975a)  has  shown,  either  approach 
would  be  of  sufficient  accuracy  with  appropriate  selection 
of  AB  and  AZ. 
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The  SSFFT  appears  to  be  juite  satisfactory 
numerically,  but  it  has  to  simplify  the  actual  oceanic  vari¬ 
ability  to  allow  realistic  computation.  It  apparently  will 
handle  source  and  receiver  depths  and  horizontal  separation 
distances  well  as  noted  by  its  degrees  of  accuracy  in  depth 
and  range.  The  smoothing  of  the  sound  speed  gradients  and 
bottom  effects  will  have  to  be  considered. 

The  intent  of  this  description  has  been  to  pro¬ 
vide  a  base  from  which  comparisons  can  oe  made  with  other 
numerical  solutions  and  modifications  to  this  algorithm. 
The  effect  of  the  source  model  on  the  MT  can  be  evaluated  by 
these  comparisons.  The  primary  elements  of  the  WT  are  the 
real  and  imaginary  parts  of  the  complex  pressure  amplitude. 
Thus  the  algorithmic  differences  with  respect  to  computation 
of  the  pressure  field  will  be  the  primary  basis  for  the  com¬ 
parison  of  effects  on  the  WT. 

McDaniel  (1975a)  compared  the  basic  Tappert  and 
Hardin  SSFFT,  which  is  used  by  Brccic,  to  normal  mode  theory. 
She  found  that  errors  arise  in  phase  and  group  speeds, 
although  one  mode  can  be  propagated  with  the  correct  ampli¬ 
tude  and  group  and  phase  speeds.  McDaniel  describes  three 
sources  of  error:  (1)  approximating  the  wave  equation  by  the 
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SSFFT  as  not<?d  above,  (2)  limiting  the  range  steps  and  per¬ 
missible  sound  speed  gradients,  and  (3)  truncating  the  field 
at  a  finite  depth.  The  result  of  the  SSFFT  will  be  a  shift 
in  the  modal  interference  at  long  ranges  if  many  modes  are 
propagating  (HcDaniel ,1975a)  .  The  error  due  to  range  step 
limitations  is  third  order  la  the  range  step  increment 
(McDaniel,  1975a)  ,  thus  a  small  range  stap  is  required  to 
minimize  this  error.  The  third  arror  is  ralated  to  inclusion 
of  a  highly  absorbing  bottom  which  causes  preferential 
attenuation  of  higher-order  modes;  the  dominant  contribution 
will  thus  result  from  the  lower  modes  (McDaniel, 1975a) . 

Brock  et  al  (1977)  have  described  a  procedure 
to  improve  this  condition.  It  maps  the  index  of  refraction 
and  depth  such  that  the  phase  speeds  and  turning  points  of 
the  normal  modes  are  preserved.  This  involves  construction 
of  a  "pseudoproblem”  such  that  the  parabolic  phase  speeds 
will  be  equal  to  the  elliptic  phase  spaads  of  the  corre¬ 
sponding  modes  in  the  original  problem  (Brock  et  al  ,  1977). 
This  is  given  as 


(n,Z)  (n,Z) 
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It  is  found  to  greatly  improve  the  SSFFT  agreement  with  the 


elliptic  solution.  This  correction  technique  has  one  basic 
constraint,  which  is  that  iaospeed  regions  are  not  allowed 
in  the  sound  speed  profile  because  of  poor  eigenfunction 
mapping  (Brock  e^  §4.  1977) . 

4.  Subsequent  Modifications 

a.  Introduction 

There  have  been  two  significant  modifications 
since  the  original  algorithm  was  published  by  Brock  (1978). 
These  are  source  function  and  volume  attenuation  alterations 
(Tatro,1977)  and  the  incorporation  of  a  bottom  and  variable 
range  step  (Stieglitz  et  al  ,  197  9).  The  Tatro  modification 
was  developed  in  1977  but  was  added  to  the  Brock  algorithm 
after  1978.  Discussion  of  rhase  two  aspects  will  complete 
the  description  of  the  PE  algorithm  installed  at  UPS  and 
utilized  in  this  study. 

b.  Tatro  Modifications 

The  Gaussian  source  originally  utilized  by  Brock 
(1978)  was  considered  to  be  inefficient  in  terms  of  program 
size  limitation  and  an  upper  frequency  limit  which  is  rela¬ 
tively  low  (Tatro,  1977)  .  A  new  source  function  was  devel¬ 
oped.  This  was  intended  to  solve  the  problem  of  large 


vertical  wavenuabers  from  a  Satis sian  source  that  were  previ¬ 
ously  eliminated  with  range  (ratro,  1977)  and  led  to  an 
incorrect  depiction  of  the  field. 

This  source  is  basically  a  low  pass  filter  in 
^he  vertical  wavenuater  domain  which  is  centered  at  the 
source  depth  (Tatro, 1 977) .  The  initial  vertical  wavenumber 
field  will  thus  be  a  constant  up  to  a  prescribed  value  and 
zero  for  higher  wavenumbers.  The  use  of  the  filter  is 
claimed  to  yield  improvement  over  the  Saussian  source  by 
minimizing  the  aliasing  which  might  occur  upon  transforming 
to  the  depth  domain  (Tatro,1977)  .  The  benefits  derived  from 
the  use  of  this  new  source  are  reduced  computation  time  and 
storage  and  the  availability  of  higher  frequencies 
(Tatro,1  977)  . 

He  also  modified  the  Broch  algorithm  by  intro¬ 
ducing  volume  absorption  as  calculated  from  the  equation  of 
Thorp  (1967).  The  decision  can  be  made  whether  or  not  to 
include  volume  absorption. 

c.  Stieglitz  gi  ^  Modifications 

Brock  (1978)  had  considered  the  PE  most  useful 
for  long  range  propagation  at  low  frequencies  and  small 
grazing  angles.  The  bottom  was  modeled  as  fully  absorbing  sc 
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that  tha  field  would  vanish  at  the  Baxiaum  depth  of  the 
transforB.  This  depth  was  obtained  by  extending  the  water 
coluan  depth  by  one- fourth  to  obtain  the  transfora  depth. 
The  index  of  refraction  in  this  region  was 


n2  =  r\2 
D 


lae 


where  a  and  b  are  empirical  constants.  The  energy  in  the 
bottom  was  attenuated  such  that  any  additional  modes  which 
resulted  from  transform  truncation  were  minimized.  This 
empirical  expression  was  selected  to  coincide  with  a  compa¬ 
rable  normal  mode  solution.  This  rather  crude  model  of  the 
bottom  appeared  to  be  acceptable  for  sitaations  wherein  the 
bottom  had  a  relatively  small  effect. 

Stieglitz  et  ^  (197  9)  considered  it  necessary 
-0  attempt  to  develop  a  more  realistic  ocean  bottom  since 
acoustic  propagation  was  being  utilized  in  more  varied  and 
extreme  situations  where  bottom  interaction  was  important. 
They  developed  two  additional  options  for  handling  the  bot¬ 
tom:  (1)  a  partially  absorbing  bottom  with  specification  of 
reflective  loss  versus  grazing  angle  of  the  equivalent  ray 
and  (2)  the  direct  insertion  of  sediment  sound  speed  and 
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attenuation  profiles.  The  algorithm  derives  the  sediment 
profiles  from  the  grazing  angle  loss  values  in  the  first 
option.  The  supplied  profiles  are  used  directly  in  the  sec¬ 
ond  option  (Stieglitz  ^  ^  ,197  9).  A  new  complex  modified 
index  of  refraction 

M(R,2)  =  {n2(R,z)  -  1}  +  i2b(Z)/K 

0 

is  computed  if  the  bottom  attenuation  is  to  be  calculated. 

One  additional  change,  the  calculation  of  the 
range  step,  has  also  been  added  by  Stieglitz  et  al  (1979) . 
It  consists  of  an  initial  search  of  the  vertical  wavenumber 
domain  for  the  maximum  pressure  component.  This  is  then  fol¬ 
lowed  by  a  search  for  the  first  component  which  is  50  dB 
below  the  maximum  component.  The  wavenumber  of  this 

component 

k  =  K  sin9 
z  o  5  0 

is  stated  to  require  a  range  step  of 

ar  =  X  /( 1  -  cose  )  , 

0  5  0 


4U 


where 


X  =  2ii/K  . 
0  0 


The  range  step  was  fixed  in  the  current  research  for  the 
benefit  of  the  spatial  transform  in  the  WI,  thus  the  above 
description  was  not  utilized.  It  has  been  included  to  com¬ 
plete  the  description  of  the  computer  coia  which  apparently 
comprises  the  NFS  model,  PEMOOEL. 

5.  SSFFT  Evaluations 

The  SSFFT  has  been  investigated  on  several  occasions 
to  examine  specific  items  which  might  be  modified  to  yield 
better  agreement  with  actual  propagarioa.  Most  of  these 
studies  deal  with  the  algorithmic  description  by  Tappert  and 
Hardin  (1973).  Examination  of  these  studies  will  note  fur¬ 
ther  possible  limitations  on  accuracy  in  the  SSFFT  which  can 
be  extended  to  the  WT . 

Fitzgerald  (1  975)  examined  the  accuracy  of  the  SSFFT 
as  a  function  of  range  and  frequency.  It  was  shown  that  the 
SSFFT  demonstrated  greater  accuracy  at  longer  ranges,  the 
lower  the  frequency.  The  SSFFT  was  considered  a  good 
approximation  in  the  examples  given  for  100  HZ  and  10  HZ  to 
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about  111  and  15000  km,  respectively  (Fitzgerald,  1975). 
The  current  study  will  ezaoine  transmission  at  50  and  100  HZ 
at  ranges  to  100  nautical  eiles.  These  ranges  will  exceed 
the  limits  of  Fitzgerald  in  some  cases. 

McDaniel  (197  5b)  investigated  splitting  the  total 
field  into  transmitted  and  reflected  fields  in  deriving  the 
SSFPT.  It  was  necessary  to  neglect  the  reflected  field  for 
the  matrix  associated  with  the  Tappert-Hardin  SSFFT,  since 
the  two  fields  did  not  decouple  where  the  wavenumber  was 
independent  of  range  (McDaniel,  1975b|  .  Thus  an  exact 
expression  could  not  be  determined  for  the  transmitted  field 
without  the  reflected  field.  (The  reflected  field  is  neg¬ 
lected  in  the  current  research  since  it  uses  the  SSFFT.) 
This  will  lead  to  phase  speed  and  group  speed  errors  as 
McDaniel  (1975a)  pointed  out  earlier,  as  discussed  in  Sec¬ 
tion  III-B-4.  These  errors  were  minimized  by  Brock  et  al 
(1977)  by  a  mapping  technique  for  both  the  index  of  refrac¬ 
tion  and  the  depth.  Even  so,  the  implication  of  rhis  possi¬ 
ble  error  cn  the  HT  is  that  the  interference  ratrerns  may 
not  be  exact,  which  will  lead  to  inexact  or  biased  wavenum¬ 
ber  nulls.  McDaniel  (1975b)  proposed  another  matrix  which 
achieved  the  desired  decoupling  and  thus  allowed  expressions 
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foe  the  transmitted  and  reflected  fields.  This  alternative 
was  not  investigated  in  this  initial  research.  It  may  be  a 
way  to  resolve  non-exact  results  and  could  be  studied  later 
to  isolate  quantitatively  the  biases  related  to  various 
input  parameters  for  the  tiT. 

Palmer  (1976)  related  possible  inaccuracies  in  the 
SSFFT  to  approximations  involving  relatively  small  terms  in 
the  exponential  approximations  for  the  operators  A  and  B  and 
also  to  the  use  of  a  single  reference  wavenumber.  These 
appreximations  will  net  be  investigated  specifically  in  this 
research  but  they  may  be  related  to  potential  disagreement 
between  the  SSFFT  and  the  theoretical  Lloyd's  Mirror  effect. 
The  association  of  errors  due  to  a  single  wavenumber  can  be 
recognized  from  the  wavenumber  spectra,  which  will  be  pre¬ 
sented  later  in  section  IV-A.  They  will  show  not  only  an 
energy  concentration  at  the  reference  wavenumber  but  signif¬ 
icant  energy  distributed,  as  well,  over  a  range  of  values. 
Palmer  (1976)  emphasized  that  the  physical  properties  of  the 
ocean  must  be  considered  before  terms  can  be  nealected,  such 
as  the  spreading  of  energy  over  several  wavenumbers, 

DeSanto  (1977)  studied  the  relationship  between  the 
acoustic  pressure  and  the  envelope  function  (called  velocity 
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potential  by  DeSanto),  which  is  directly  relevant  to  the 
current  research  where  pressure  is  calculated  froa  the  aod- 
eled  envelope  function.  DeSanto  arrives  at 


iJi(R,Z)  =  ^"P(t,2)  Q(R,t,Z)  exp  ( iK^/Zt)  (R2  +  t2) 


where  Q  satisfies 


«RR  *  ♦(iK„R)/t)Qp  *  + 


2iK^Q^  *  2(!^^(ln<p(t,Z)>) 

32 


K^2(K2(t,Z)  -  K^CR^Z)  }Q 


and 

K2(R,Z)  =  K(R,Z)  -  K^(Z)e. 


The  PE  for  the  envelope  function  is  obtained  if  this  func¬ 
tion  is  evaluated  by  stationary  phase;  thus,  it  appears  that 
a  acre  accurate  result  could  be  obtained  by  better  evalua¬ 
tion  of  the  integral.  DeSanto  et  (1978)  developed  a  cor¬ 
rection  for  the  SSPFT,  known  as  the  Corrected  Parabolic 
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Approxiaation  (CPA),  that  handled  phase  inaccaracies  better 
than  the  SSFFT.  Further  study  might  address  the  response  of 
the  HT  to  the  CPA  with  respect  to  improved  phase  accuracy  as 
a  function  of  range. 

6.  laaiigaBtaiai  of  thf  pe 

This  final  section  in  the  analysis  of  the  PE  will 
consider  the  expected  environmental  sensitivity  of  the  PE, 
and,  by  extension,  the  iT.  The  important  relationship  of  the 
PE  to  the  iT  is  the  input  of  the  complex  pressure  amplitude 
and  ultimately  the  resultant  wavenumber  distribution.  As 
discussed  earlier  in  Section  11- B,  the  wavenumber  could  be 
associated  to  the  source  angle  and,  by  extension,  to  the  ray 
type  or  family  (Lauer,  1979).  Consequently,  oceanic  parame¬ 
ters  which  modify  the  transmission  of  cartain  rays  can  be 
related  to  the  resulting  wavenumber  distribution. 

The  sound  speed  profile  will  favor  certain  ray  paths 
at  different  depths  in  the  water  column.  There  should  be  a 
corresponding  shift  in  ray  paths  and  hence  wavenumbers  as 
the  source  and  receiver  depth  geometry  is  changed.  The 
direct  and  surface-reflected  rays  will  travel  at  different 
speeds  depending  upon  the  gradients  along  the  profile  and 
thus  lead  to  different  spatial  interference  patterns.  This 
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can  leal  to  wavenuaber  shifts  aa d  changes  in  the  null  spac* 
ing.  Profiles  used  in  the  current  study  will  include  real 
exaaples  and  isospeed  cases. 

Another  important  aspect  is  the  transmission  geom¬ 
etry.  Generally  the  model  will  be  run  at  considerable  range 
and  in  relatively  deep  water,  which  will  be  consistent  with 
the  applicability  of  the  PE.  A  series  of  source  depths  with 
a  single  fixed  receiver  depth  will  be  run  to  generate  the 
wavenumber  null  spacing  values  that  will  formulate  a  deter¬ 
mination  curve  for  each  case  of  geometrioal  variation.  The 
distribution  of  these  resulting  curves  can  be  examined  to 
evaluate  the  effects  cf  transmission  geometry  on  the  WT. 

Two  frequencies  will  be  examined  in  a  similar  man¬ 
ner.  This  would  be  one  of  the  most  important  parameters  to 
investigate  in  future  studies  since  it  would  have  a  direct 
bearing  on  operational  applications.  It  is  expected  that 
there  will  be  differences  in  the  resulting  wavenumber  spec¬ 
tra  depending  on  the  input  frequency.  The  relationship 
between  frequency  and  transmission  geometry  will  also  be 
studied  since  these  are  the  easiest  parameters  to  vary  dur¬ 
ing  an  experiment  and  will  lead  to  the  most  complete  analy¬ 
sis  of  model  variability. 
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The  surface  will  be  assumed  to  be  planar  and  pres¬ 
sure-release  so  that  there  will  be  complete  specular  reflec¬ 
tion.  The  minimum  source  depth  used  will  be  that  at  which  a 
definable  interference  pattern  can  be  observed  for  the  spec¬ 
ified  frequency.  The  slope  approaches  zero  near  the  surface 
on  the  source  depth  determination  curves.  This  result  may 
indicate  that  the  WT  in  its  current  form  is  probably  not 
applicable  for  sufficiently  shallow  sources. 

The  final  important  area  of  variability  will  be  the 
bottom  boundary  conditions.  There  are  several  options  now 
available  in  the  PEMODEL.  The  fully  absorbing  and  partially 
absorbing  bottoms  will  be  considered.  k  fully  absorbing 
bottom  will  be  used  in  most  cases  to  simplify  the  interfer¬ 
ence  pattern  by  elimirating  any  rays  which  interact  with  the 
bottom.  Thus  the  basic  interference  pattern  may  be  evalu¬ 
ated.  Various  bottoms  can  then  be  tested  to  study  the 
effect  of  the  bottom  on  the  rfT.  Differing  bottom  less 
curves  will  be  compared  in  a  few  cases  to  demonstrate  quali¬ 
tatively  that  the  bottom  does  have  an  effect.  The  third  bot¬ 
tom  cption  of  a  realistic  sediment  sound  speed  could  also  be 
investigated  once  the  general  effect  of  the  bottom  has  been 
determined.  It  would  not  offer  any  additional  significant 
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results  at  this  point  beyond  those  realized  from  the  fully- 
and  partially-  absorbing  bottoms. 

All  input  parameters  are  related  and  thus  the  exami¬ 
nation  of  any  one  particular  factor  must  be  considered  with 
respect  to  the  full  range  of  variation  of  the  other  factors. 
The  number  of  combinations  of  parameters  would  be  very  large 
for  a  model  such  as  the  PEHODEL.  Thus  this  research  can 
only  begin  to  explore  general  relationships.  Additional 
research  is  necessary  to  delineate  specific  parameter 
impacts  on  the  WT. 

C.  FAST  FIELD  PROGS  AH 

1 ,  General  Descr  iption 

The  Fast  Field  Program  (FFP)  was  selected  as  the 
first  model  to  compare  to  the  PEMODEL  since  it  was  the  ori¬ 
gin  of  the  WT  and  the  initial  testing  of  the  WT  by  Lauer 
(1979)  was  performed  with  the  FFP.  The  description  of  the 
FFP  which  follows  is  summarized  from  DiSapoli  (1971). 

The  basic  result  of  the  FFP  is  direct  numerical 
integration  by  the  application  of  the  FFT  to  field  theory  to 
compute  propagation  predictions  in  a  minimum  computation 
time.  The  pressure  can  be  represented  in  the  form 

?(R,Z)  =  /T"/”  G(Zj^,Z2;K)  H^^KR)  K  dK  ^ 
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where  G,  the  Green's  function,  aust  satisfy 


d^G  +  {K2(Z)-k  2}g  =  -5(Z-Z  ) 
-  r  s 

dz^ 


and  the  associated  boundary  conditions.  The  field  equation 
can  then  be  written  as 


Z,R  ) 
’  m 


=  fa  AK(  2/Tri)  ^  e 


iK  R 
0  m 


/R  ^  E 
m  rti 


LuR  AK 
e  0 


where 


=  G(Z,Zc.;K^)  K 
m  bra. 


ra 


\  imR  AK 
e  0 


In  this  form,  p  can  be  obtained  through  use  of  the  FFT  at 
each  range  increment.  Given  p  as  a  function  of  range,  the  W'T 
could  be  applied  as  usual.  DiNapoli  compared  an  example 
from  the  FFP  with  the  normal  mode  prograa  of  Bartberger  and 
Ackler  (1973)  and  found  excellent  agreement  between  the 
locations  of  the  peaks  and  the  real  parts  of  the  eigenval¬ 
ues.  The  only  significant  difference  is  that  the  FFP 
includes  more  higher  order  modes  but  does  not  _nclude  the 
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first  four  nodes.  DiNapoli  states  that  the  FFP  offers  a 
significant  reduction  in  computation  time  yet  still  provides 
reasonable  results  as  compared  to  normal  lode  theory.  This 
is  accomplished  by  expressing  the  sound  velocity  profile  as 
one  or  more  exponential  functions  of  depth.  These  exponen¬ 
tial  functions  allow  the  use  of  recurrence  relations  to 
guickly  calculate  the  input  to  the  FFT. 

The  FFP  offers  two  options  for  the  modeling  of  the 
bottom;  a  two-layered  fluid  bottom  wi'-h  specified  sound 
speeds  and  a  semi-infinite  bottom  of  one  sound  speed.  The 
propagation  effect  from  the  first  would  be  eventual  return 
of  rays  once  reflected  from  the  second  bottom  interface.  The 
semi-infinite  bottom  would  not  allow  any  return  of  sediment 
ra  ys. 

2.  Comparison  of  FFP  and  PS 

The  FFP  is  a  range  independent  model  and  the  PE  is 
range  dependent.  The  environment  was  considered  constant 
over  the  range  in  the  current  research,  and  thus  the  two 
models  can  be  compared.  The  FFP  uses  the  sound  speed  profile 
as  an  exponential  function  which  is  a  different  approach 
from  the  linear  segmentation  of  the  PE.  The  number  of  expo¬ 
nential  functions  could  be  sufficient  to  model  an  actual 
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profile  approxiaately  as  well  as  the  PE.  The  PE  bottom  mod¬ 
eling  does  net  appear  to  offer  any  significant  advantage  in 
theory  over  the  PPP,  which  allows  for  a  semi-infinite  or 
layered  bottom.  It  is  likely  from  a  theoretical  standpoint 
that  the  PE  and  FFP  inputs  to  the  KT  could  be  similar  enough 
to  allow  comparison.  The  wavenumber  spectrum  shows  a  series 
of  peaks  which  result  from  the  nearby  singularities  of  the 
Green's  function.  These  can  be  expressed  as  the  normal 

modes.  The  utility  of  the  WT  stems  from  the  spacing  of  the 
nulls  between  the  envelopes  which  contain  these  peaks. 

0.  NORMAL  MODE  THEORY 

1 .  General  Deset  IptirP?^ 

The  fundamental  importance  of  the  normal  mode  solu¬ 
tion,  and  the  reason  that  is  it  used  no  check  the  validity 
of  other  models,  is  that  it  is  an  exact  solution  to  the  wave 
equation  (Kir.sler  et  gl  ,  1982 ,p.  430|  .  The  following  descrip¬ 
tion  continues  from  Kinsler  et  ^1  (1982,  pp.  U30-432) . 

The  solution  for  a  point  source  is  given  by 

{ v2-32/c23t2)  0  =  1  , 

4  ir 

where 

^0  =  • 
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Analysis  yields 


p  =  (-jTr)ZZ„(2«)Z„(Z)H  UK  R) 

■'nnSn  on 

and  Zn  oust  satisfy 

d^Z„  +  {u2/c2(Z)  -  K  2}Z  =  0  , 

_ n  n  n  * 

dz^ 

where  Kn  is  constant  and  appropriate  boundary  conditions. 
This  solution  is  for  trapped  modes  and  does  not  include 
those  which  are  evanescent.  normal  mode  theory  can  be 
linked  to  the  WT  by  the  spatial  representation  of  the  wave- 
number 

K(Z)  =  _u) _ 

C(Z) 

The  rays  in  the  deep  sound  channel  which  have  the  same  local 
direction  of  propagation  correspond  to  the  normal  modes  by 
the  relationship 

COS0  =  , 

K(Z) 
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which  extends  to  the  normal  aode  phase  speed 


where  1  designates  the  1th  moda. 

Normal  mode  theory  lands  itself  to  many  applications 
in  underwater  acoustics  and  it  may  be  formulated  for  both 
range  dependent  and  independent  situations,  k  dependent  form 
which  is  appropriate  for  comparison  with  the  SSFFT  was 
developed  by  Kanabis  (1975).  He  developed  a  normal-mode 
model  which  allows  large  changes  in  the  sound  speed  profile, 
depth,  and  bottom  composition  with  ranga.  The  range  domain 
is  segmented  into  regions  within  each  of  which  sound  speed 
and  bottom  composition  are  functions  of  depth  but  not  range. 
The  normal  mode  solutions  for  each  region  are  matched  at  the 
interfaces. 

2.  Coj8£a£i§ofl  ol  NM  a^  £E 

It  would  be  expected  that  the  NT  would  respond  simi¬ 
larly  to  these  two  models.  This  is  based  upon  this  example 
of  a  normal  mode  program  and  supported  by  the  previously 
discussed  evaluation  studies  involving  PE  and  NM.  The  NM 
solutions  are  generally  applicable  to  long  range  propagation 
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at  lc¥  frequencies,  which  is  quite  similar  to  the  PE.  The 
NM  can  be  identified  closely  with  the  HT  in  a  physical 
sense,  because  of  the  interference  effects  in  the  NH  that 
result  froB  several  multiple  reflections  (Officer,  1958, p. 
117),  Officer  further  states  that  the  nuaber  of  reflections 
will  increase  and  the  time  intervals  between  the  incidence 
of  successive  reflections  will  decrease  as  the  range 
increases.  This  might  be  comparable  to  depth  variations 
which  will  be  shown  for  the  WT.  The  SM  yields  exact  solu¬ 
tions  for  modeling  with  planar  boundaries.  Thus  the  results 
of  a  WT  based  on  NM  would  be  limited  only  by  approximations 
which  were  intrinsic  to  the  WT.  A  comparison  of  PE-  and  NM- 
based  WT  data  should  aid  in  evaluating  the  impact  of  the 
particular  PE  algorithm, 

E.  FINITE  DIFFERENCE  METHODS 
1 .  General  De^r iptions 

One  final  group  of  methods  which  has  evolved 
recently  are  the  finite  difference  methods  of  solution  for 
the  parabolic  equations  (Lee  and  Papadalcis,  1979).  They 
seek  a  more  general  solution  which  will  be  appropriate  in 
shallow  water  or  where  bottom  interaction  will  be  important. 
They  further  state  that  the  PD  and  ODE  methods  will  be 
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superior  to  SSFPT  tecause  an  artificial  bottom  is  not 


required  for  the  benefit  of  tha  PPT.  Explicit  and  implicit 
finite  difference  schemes  are  considered  as  well  as  an  ordi¬ 
nary  differential  equation.  The  explicit  finite  difference 
(EPD)  scheme  is 

n+l  K  .  n 
u  =  e  3  u 

m  —  m  * 

3R 


where  u^  represents  4/(Rn,2m),  8n  is  the  nth  range  point,  and 
Zn  is  the  mth  depth  point.  The  implicit  finite  difference 
(IPD)  method  from  Lee  and  Botssas  (1982)  is 


e  3  R  u. 


—  n+l 


HK3/3R 


m 


=  e 


u 


m 


The  FD  scheme  is  consistent,  stable  and  convergent  in  the 
solution  of  the  PE  (Lee  and  Papadakis,  1979). 

The  ODE  method  is 


d  u_  =  a_  u_  +  b»«  (u  .  -2u  +u  ) 

—  m  mm  -  m+ 1  m  m-l  > 

dR  (AZ)2 


where  a  and  b  are  related  by 


u  =  a(K^,R,Z)u  +  b(K  ,R,Z)u, 

Is  0  0 
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This  technique  is  stable,  consistent  and  convergent  also 
which  supports  its  use.  One  additional  aspect  of  the  ODE 
approach  is  the  calculation  of  a  variable  range  step. 

2.  Coaparison  al  P£  aijd  ODE  with  gE 

One  result  of  the  FD  and  ODE  aethods  is  greater 
flexibility  in  arbitrary  bottom  and  surface  boundary  condi¬ 
tions  than  that  available  with  the  SSFFT.  Lee  and  Papadakis 
(1979)  state  that  while  the  FD  and  ODE  methods  are  both 
superior  to  the  SSFFT,  the  best  appears  to  be  the  ODE,  fol¬ 
lowed  by  the  EFD,  which  is  limited  by  step  size.  Thus  the  FD 
and  ODE  methods  appear  to  be  significant  improvements  over 
the  SSFFT.  This  result  should  extend  to  a  HT  based  upon  FD 
or  ODE  input. 

Wore  recently  Lee  et  al  (1981)  have  re-examined  the 
IFD  and  found  it  to  be  as  equally  efficient  as  the  DDE.  The 
general  advantage  of  the  IFD  over  the  SSFFT  is  that  the 
problem  is  solved  within  the  water  column  since  the  bottom 
effects  can  be  modeled  by  straight  line  segments.  This  is 
similar  to  the  ODE,  but  the  IFD  does  not  require  the  amount 
of  storaqe  necessary  for  *he  current  ODE  (Lee  et  3.I,  1981). 
McDaniel  and  Lee  (198  2)  have  extended  the  IFD  to  the  treat¬ 
ment  of  vertical  density  discontinuities.  They  found  that 
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while  without  this  eore  realistic  treataent  of  the  water- 
bottOB  interface  the  IFD  approriaates  ths  SSFPT  and  differs 
significantly  fro«  noraal  aode,  the  IFD  corresponds 
extreaely  well  with  the  normal  mode  calculations  with  the 
interface.  In  conclusion^  the  numerical  techniques  related 
to  FD  and  ODE  appear  to  be  viable  alternatives  to  the  SSFFT. 
They  will  permit  a  more  accurate  representation  of  underwa¬ 
ter  acoustic  transmission,  which  will  be  directly  relatable 


11.  ANALYSIS  OP  NT 


A.  GENERAL  DESCRIPTION 

Twenty  scenarios  were  selacted  to  eranine  the  general 
aspects  of  variations  in  several  oceanic  and  geonetric 
paraaeters.  A  set  of  input  parameters  and  a  sound  speed 
profile  were  constant  in  each  scenario  while  the  source 
depth  was  varied.  Five  to  ten  source  depths  were  used  for 
each  situation.  The  vertical  and  horizontal  axes  of  the 
wavenumber  spectra  correspond  to  the  spectral  intensity  and 
■^he  scaled  wavenumber  beta,  respectively.  The  intensify 
axis  was  normalized  to  unity  for  ease  in  plotting.  The  null 
spacing  was  determined  for  each  spectrum.  These  null  spac¬ 
ing  values  were  plotted  versus  source  depths  to  yield  the 
determination  curve  for  each  scenario. 

Table  I  lists  the  paramerric  specifications  for  each 
scenario.  Basically  they  can  be  divided  into  isospeed  or 
sound  speed  profile  cases.  An  isospeed  profile  of  4890  ft/s 
from  the  surface  to  the  bottom  was  seieC-ed  *0  evaluate  the 
transmission  geometry  and  the  general  consistency  of  experi¬ 
mental  results,  assuming  a  fully-absorbing  bottom.  The  <JT 


is  based  on  the  arrivals  of  the  direct  and  surface-reflected 
rays;  thus,  the  Isospeed  cases  will  be  the  easiest  to  con- 
pare  with  theory.  There  will  be  no  sound  speed  differences 
as  a  function  of  depth  which  could  aodify  the  arrival  tines. 
Additionally,  several  realistic  ocean  profiles  were  studied 
to  allow  a  qualitative  coaparison  with  theory  and  explore 
spatial  and  teeporal  variability.  It  is  recognized  that  the 
theoretical  Lloyd's  Mirror  effect  applies  only  to  isospeed 
profiles.  Besults  from  real  profiles  were  compared  with  this 
theory  to  see  if  they  could  approximate  it.  The  spatial  and 
tempor2d  scales  of  variability  used  in  this  research  were 
approximately  200  miles  and  20  days. 

Bore  analyses  were  performed  with  the  isospee’d  cases 
since  this  was  an  initial  investigation.  It  still  remains 
to  consider  other  simple  cases,  such  as  Linear  negative  and 
positive  gradients,  surface  ducts  and  deep  sound  channels. 
Also  realistic  ocean  studies  that  would  involve  greater  tem¬ 
poral  and  spatial  variations  would  be  beneficial.  Time  was 
limited  during  rhis  research,  thus  efforts  were  divided 
between  model  validation  and  oceanic  sensitivity.  It  is 
believed,  however,  that  the  comparisons  which  will  be 
presented  are  sufficient  for  a  preliminary  investigation 
into  the  sensitivity  of  the  WT. 


T&BLS  I 

Scenario  Specifications 
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NuBber  of 

Set 

HD:  ft 

3  :nB 
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Z:ft 

SD  Runs 

10 

300 

50 

100 

10000 

7 

11 

300 

50 

50 

10000 

6 

12 

300 

100 

50 

10000 

6 

13 

300 

50 

100 

5000 

6 

14 

300 

100 

100 

5000 

5 

15 

300 

50 

50 

5000 

6 

16 

300 

100 

50 

10000 

6 

17 

300 

100 

100 

10000 

5 

18 

800 

50 

50 

10000 

5 

19 

10000 

50 

50 

10000 

5 

HI 

300 

25 

100 

10000 

5 

32 

300 

50 

100 

10000 

5 

33 

300 

75 

100 

10000 

5 

REALISTIC  PROFILE 

Bun 

Rev  Depth 
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Frsq 

Rater  Depth 

Sre  Depth 
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SD  (number! 

AO 

300 

100 

100 

9000 

8 

A1 

1000 

100 

100 

9000 

5 

A2 

9000 

100 

100 

9000 

5 

A3 

300 
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100 

9000 

5 

BO 

300 

100 

100 

12000 

8 

CO 

300 

100 
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9000 

8 

00 

300 

100 

100 

12000 

8 
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The  coaplex  pressure  input  data  were  not  modified  so 
that  the  entire  wavenumber  spestrue  would  be  observed  at  the 
receiver.  These  effects  would  need  to  be  evaluated  if  the 
HT  was  to  be  compared  to  received  signals.  The  surface  and 
bottom  were  assumed  flat.  Boundary  variations  would  have  to 
be  considered  for  WT  use  under  realistic  ocean  conditions. 
A  nominal  beam  width  of  thirty  degrees  was  used  and  in 
almost  all  cases  the  range  step  was  3.1  mile. 

frequencies  of  50  and  100  Hz,  water  column  depths  of 
5000  and  10,000  feet  and  ranges  of  50  and  100  nautical  miles 
were  the  basic  comparisons  for  most  isosoeed  scenarios.  A 
receiver  depth  of  300  feet  and  a  fully  absorbing  bottom  were 
utilized,  unless  otherwise  noted.  Additional  comparisons 
were  performed,  such  as  range  variations  to  25  and  75  miles 
and  receiver  depth  variations  of  000  and  10,000  feet.  A 
range  of  100  nautical  miles  and  a  frequency  of  100  Hz  were 
used  for  the  actual  cases  unless  otherwise  noted.  The  bottom 
depths  were  9000  feet  for  profiles  A  and  C  and  12,000  feet 
for  profiles  5  and  0. 

The  realistic  profiles  were  obtained  from  the  Acoustic 
Storm  Transfer  and  Response  Experiment  (ASTREX)  which  was 
conducted  by  NFS  in  the  northeast  Pacific  Ocean  in  November 
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and  OecsBbec,  1980  (Ounlap,1982)  .  Ihe  ASIREX  study  vas  con* 
ducted  to  evaluate  the  effect  of  oceanic  rssponse  to  atmos¬ 
pheric  storms  or.  acoustic  propagation.  The  profiles  were 
selected  for  this  research  to  represent  only  spatial  and 
temporal  variability  without  respect  to  any  experiment  oper¬ 
ations.  The  profiles  were  obtained  from  expendable  bathy¬ 
thermographs  (AXBT)  which  were  dropped  by  0.  S.  Navy  P-3 
aircraft.  Fig.  4  shows  the  relative  locations  and  times  for 
these  profiles.  The  digitized  profiles  were  analyzed  by  the 
Oniversity  cf  Hawaii  to  remove  spikes  and  false  values. 
These  profiles  were  combined  with  climatological  salinity 
profiles  at  the  Fleet  Numerical  Oceanography  Center  in  Mon¬ 
terey.  Table  II  gives  the  descriptive  data  for  each 
profile.  The  bottom  loss  curves  which  were  used  are  similar 
to  those  found  in  Orick  (1979)  and  were  arbitrarily  devel¬ 
oped  to  indicate  higher  loss  for  profiles  3  and  D.  This  is 
related  to  bottom  composition;  bottom  types  for  profiles  A/C 
and  B/D  are  calcareous  sand  and  clay,  respectively 
(NAVOCBANO,1 978)  .  Thus  higher  loss  is  expected  for  B/D. 

A  general  feature  of  the  spectra  for  almost  all  test 
cases  is  a  distinct  0-shape  with  peak  values  at  the  minimum 
and  maximum  wavenumbers.  Figs.  5  through  10  show  this  phe- 
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TABL2  II 

Location  Information  for  Realistic  Profiles 


Profile 

Lat 

Long 

Date 

Time (Z) 

NLD(m) 

A 

U1-2aN 

127-23i 

16  NOV 

1928:04 

32 

B 

45-47N 

133-05H 

16  NOV 

2121:25 

61 

C 

41-23N 

128- 22H 

6  DEC 

1843: 37 

58 

D 

45-46N 

133-09H 

S  DEC 

2033:45 
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r.omenon  for  isospeed  cases  and  Figs.  11  through  16  for  a 
realistic  sound  speed  profile.  The  abbreviations  on  these 
figures  are  source  depth  (SO),  receiver  depth  (RD)  ,  and 
water  column  depth  (Z)  in  feet;  range  (R)  in  nautical  miles; 
and  frequency  (F)  in  Hz.  This  key  will  also  apply  to  the 
source  depth  determination  curves  which  will  be  presented 
later.  The  righ-^.  and  left  intensity  maxiaa  appear  no  corre¬ 
spond  to  beam  eleveation  angles  of  0  and  30  deg,  respec- 
•^ively,  which  could  be  considered  to  approximate  the  direct 
and  surface-reflected  waves.  A  30  deg  beam  angle  was  used; 
thus,  there  may  be  a  relationship  with  the  left  maximum. 
This  beam  angle  is  also  near  the  maximum  angle  which  Tappert 
(1977)  considered  to  be  appropriate.  Then  the  l“ft  maximum 
might  be  a  result  of  algorithmic  inaccuracy.  Different  beam 
angles  could  be  tested  to  determine  if  this  would  explain 
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the  unexpected  D-shape  curve.  This  apparently  anoaalous 
spectral  pattern  might  be  explained  by  a  possible  transform 
assumption  of  symmetry  about  zero  and  the  inclusion  of  neg¬ 
ative  transform  values.  It  is  also  possible  that  the  left 
wavenumber  peak  may  be  the  result  of  aliasing  of  the  high 
wavenumber  values.  OiHapoli  (1971)  suggested  that  aliasing 
could  occur  and  thus  he  utilized  only  one-half  of  the  range 
field.  This  has  not  been  tested  in  the  current  research; 
however,  examination  of  different  range  segments  might 
resolve  this  apparent  problem.  A  remote  possibility  is  that 
the  direct  and  surface-reflected  rays  each  exhibit  a  Bessel 
function  display.  This  might  be  depicted  by  a  peak  at  the 
reference  wavenumber  for  each  ray  with  a  set  of  interfering 
envelopes  between  the  peaks.  A  mathematical  explanation  for 
this  possibility  has  not  yet  been  explored.  No  spikes  are 
seen  in  very  shallow  water  as  depicted  in  Fig.  17,  but  also 
no  distinct  null  pattern  is  observed,  thus  near  surface 
source  depths  have  been  neglected. 

These  figures  were  also  included  to  show  a  series  of 
source  depth  runs  and  the  resulting  decrease  in  bet?  spacing 
as  a  function  of  source  depth.  The  null  spacing  was  manu¬ 
ally  measured  for  each  run  as  the  distanca  between  two  adja- 
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NQimiZCO  SPCCTRRL  INTENSITY 
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NQRfmLIZCQ  SPCCTRflL  INTENSITY 
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cent  intensity  minima,  such  as  shown  in  Figs,  5  through  16. 
The  nulls  selected  for  measurement  were  rhose  in  the  middle 
of  the  spectrum,  which  led  to  a  consistent  technique  for 
comparison  and  avoided  the  rapidly  decreasing  null  spacing 
near  the  left  intensity  maxima.  The  distance  was  measured 
to  the  nearest  1/64  of  an  inch  and  equated  to  the  respective 
null  spacing.  Thus,  there  is  some  inheren*  measurement 
error.  The  exact  null  locations  were  difficult  tc  determine 
in  the  shallow  isospeed  cases  because  of  ger.'^  le  spectrum 
curvature  and  uncertainty  in  null  location.  The  actual 
profiles  were  very  complex  and  again  the  nulls  were  diffi¬ 
cult  to  locate.  One  value  for  each  spectrum  plot  was  used 
'TO  generate  a  point  on  the  determination  curve,  such  as  Fig. 
18  for  a  series  of  rcns.  The  variation  of  the  data  poin-.s 
from  the  theoretical  Lloyd’s  Mirror  curve  will  be  the  normal 
distance  between  the  point  and  the  curve.  Measurement  inac¬ 
curacy  will  decrease  as  source  depth  increases.  This  will 
be  significant  because  at  shallow  depth  there  is  some  varia¬ 
tion  noted  from  sound  speed  differences.  The  nulls  are 
well-defined  at  great  depths,  but  there  is  lit-le  varia'^ion 
already,  thus  the  measurement  errors  were  net  considered 
significant.  It  is  possible  that  a  numerical  algorithm  which 
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could  quantitatively  identify  the  null  locations  could  mini¬ 
mize  this  error. 

This  expectation  is  inferred  from  the  wavenumber  plots 
by  Lauer  (1979)  and  his  examination  of  range  determination. 
He  found  that  range  determination  resolution  increased  as 
the  distance  between  source  and  receiver  decreased.  It  is 
possible  that  a  similar  relationship  might  exist  for  depth 
determination.  Thus  greater  resolution  aay  be  realized  at 
shallow  source  depth.  Even  at  greater  source  depths  the  res¬ 
olution  may  be  acceptable  for  certain  applications  where  an 
approximate  depth  is  acceptable. 

The  prediction  curves  which  result  from  the  individual 
HT  spectra  will  be  discussed  in  following  sections.  Table 
III  is  a  listing  of  the  prediction  curve  comparisons  which 
were  conducted.  in  the  figures  only  one  symbol  may  appear 
at  a  source  depth  if  both  curves  had  the  same  null  spacing 
value. 
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TABLE  III 

Coaparison  Specifications 


ISOSPEEO 

PROFILE 

Figure 

Run  Sets 

Constant  Parameters 

Varying  Parameter 

19 

11,15 

R:  50,F:50,  RO:  300 

2:10000,5000  ft 

20 

10,13 

R:50,F:100,RD:300 

2:10000,5000  ft 

21 

12,16 

R:  100,F:50,RD  :300 

2:10000,5000  ft 

22 

17,14 

R:  100,  F:100,  RD:300 

2:10000,5000  ft 

23 

15,16 

F: 50,2:5000, RD; 300 

R:50,100  nm 

24 

11,12 

F:  50, 2:10000, RD:300 

R:50, 100  nm 

25 

R1,R2,R3 

F:  100,2:10000  ,  RD:  300 

R:25,50,75  nm 

26 

12,17 

R :  100 , Z : 10000 , RD : 300 

P:100,50  h2 

27 

10,11 

R:  50,2:  10000,  RD:300 

F:100,50  hz 

28 

11^18,19 

R;  50,  F:50,2:  1  0000 

RD; 300,800, 10000  ft 

REALISTIC  PROFILE 

Figure 

Run  Sets 

Constant  Parameters 

Varying  Parameter 

29 

AO,  BO 

R;  100,F:100,RD:300 

Space: 200  nm 

30 

CO, DO 

R:  100,F:100,RD:300 

Space:200  na 

31 

AO, CO 

R:  100,  F:100,rD:300 

Tim9:20  days 

32 

BO, DO 

R:  100,F:100,rD:300 

riine:20  days 

33 

AO,  A3 

F:  100,  16N0V,  RD:300 

8:100,50  nm 

34 

A0,A1,A2 

F:  100, 16NOV,  R;ioQ 

RD;  300, 1000, 9000  ft 

COMPOSITE 

Figure 

Hun  Sets 

Constant  Parameters 

Varying  Parameter 

35  (All  Isospeed  and  Actual  Profile  Run  Sets) 
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B.  ISOSPEED  C&SBS 


1-  Salaia 

Pig.  19  is  the  coaparison  for  a  frequency  of  50  Hz 
and  a  range  of  50  na  at  vater  column  depths  of  5000  and 
10,000  feet.  It  was  not  expected  that  there  would  be  any 
variation  if  the  water  cclumn  depth  was  modified,  since  all 
bottca  interacting  rays  should  be  fully  absorbed.  The 
curves  showed  excellent  agreement  at  source  depths  of  800 
feet  and  greater,  but  there  was  some  variation  noted  for  500 
feet  and  shallower.  Some  variation  may  be  related  to  meas¬ 
urement  error,  as  noted  above.  Maximum  variation  from 
*heoty  is  observed  for  both  depths  in  the  middle  depth 
range,  approximately  500  feet,  which  is  also  the  region  of 
greatest  curvature  in  the  source  depth  determination  curve. 

Fig.  20  is  the  same  except  the  frequency  has  bean 
increased  to  100  Hz.  There  is  better  agreement  from  200  to 
3000  feet;  however,  the  significance  of  the  frequency  varia¬ 
tion  has  not  been  evaluated.  Fig.  21  is  for  a  frequency  of 
50  Hz  and  a  range  of  100  nm,  and  Fig.  22  is  for  100  Hz  and 
100  nm.  Some  variaticns  are  noted;  however,  they  may  not  be 
significant,  which  would  agraa  with  the  supposition  that 
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)  SOURCE  DEPTH  lET) 

1000.0  tSOO.O  2000.0  zst 


SOURCE  DEPTH  OETERM I NATION  CURVE 
ISQSPEEO  WATER  DEPTH  VARIATION 


acm  NUi.  sPHCiw  (i/rri 


Figure  19 


sospeed  Profile  SD  Detarmination  Curve:  Z 
ariition  •  1. 


water  coluon  depth  differences  will  not  be  significant  for 
an  isospeed  profile. 

2.  Sssae 

Pig.  23  is  the  cosparison  for  a  frequency  of  50  Hz 
and  a  water  colunn  depth  of  5000  feet  for  50  and  100  na. 
There  is  essentially  no  differsnce  between  the  curves.  The 
experimental  curves  agree  well  from  800  to  3000  feet,  but 
once  again  the  shallow  variations  may  result  from  inaccura¬ 
cies  in  evaluating  beta  from  the  spectrum  plots.  Pig.  24  is 
a  similar  comparison,  but  the  water  column  depth  has  been 
increased  to  10,000  feet.  Fig.  25  increases  the  frequency 
to  100  Hz  for  ranges  of  25,  50  and  75  miles.  Again  all 

curves  show  very  good  agreement  with  theory,  which  indicates 
that  range  variations  within  100  nm  of  the  source  will  prob¬ 
ably  be  insignificant.  Ranges  beyond  100  nm  were  not  tested 
since  Fitzgerald  (1975)  suggested  that  the  SSFFT  was  not  a 
good  approximation  for  100  HZ  beyond  111  Icm.  This  could  be 
tasted  later. 


3.  Frequency  Yar  igit^on 


Fig. 

26  is 

the 

comparison 

for 

a  range  of 

50  nm 

and  a 

water  column 

depth 

of 

5000  feet 

for 

50  and  100 

Hz. 

Both 

curves  agree 

well 

with 

theory  below 

300  feet. 

but 

at  500 
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SOURCE  DEPTH  DETERMINATION  CURVE 
ISOSPEED  WATER  DEPTH  VARIATION 
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Figure  20.  Iso^peed  Profile  SD  Detsrmiaation  Curve 
Variation  -  2. 
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SOURCE  DEPTH  DETERMINATION  CURVE 
ISOSPEED  WATER  DEPTH  VARIATION 


ecm  NULL  SPACING  (1/PT) 


Figure  21.  Iso^peed  Profile  SD  Determination  Curve 
Variation  -  3. 
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SOURCE  DEPTH  DETERMINATION  CURVE 
. ISOSPEED  WATER  DEPTH  VARIATION 


22.  ^|0$ge|d^ Profile  SD  De'tsrmination  Curve; 


SOURCE  DEPTH  OETERMINRTION  CURVE 
ISOSPEED  RANGE  VARIATION 
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Iso$Deed  Profile  SD  Determination  Curve 
Variation  -  1. 


Figure  23. 
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SOURCE  DEPTH  OETERMINRTION  CURVE 
ISOSPEED  RANGE  VARIATION 
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Figure  24.  Isosoeed  Profile  SD  Detarminat ion  Curve;  R 
Variition  2. 
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SOURCE  DEPTH  DETERMINATION  CURVE 
ISOSPEED  RANGE  VARIATION 
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Fiqure  25.  Isospeed  Profile  SD  Determination  Curve:  R 
Variation  -  3. 
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feet  the  higher  frequency  appears  to  agree  better.  Fig.  27 
extends  the  water  colusn  depth  to  10,000  feet  and  there  is 
apparently  betc-vc  agreement  with  the  theoretical  curve. 
Overall  there  appears  to  be  better  agreeient  at  100  Hz  than 
at  50  Hz,  but  this  result  is  not  certain  with  the  limited 
number  of  cases  run. 

4.  Receiver  Depth  Variation 

Fig.  28  compares  receiver  depths  of  3C0,  8  00  and 
10,000  feet  for  a  range  of  50  nm,  water  column  depth  of 
10,000  feet,  and  frequency  of  50  Hz.  Again  all  curves  agree 
reasonably  well.  There  is  perhaps  slightly  better  agreement 
for  the  shallow  receiver,  but  this  supposition  cannot  be 
verified  from  these  limited  data. 

5.  An alvsis  o£  Isospeed  Cases 

There  do  not  appear  to  be  any  significant  differ¬ 
ences  for  water  column  depth,  range,  frequency  or  receiver 
depth  over  the  limits  of  variations  considered.  This  was 
expected  since  the  theory  is  supposed  to  be  only  a  function 
of  source  depth  for  an  isospeed  profile.  Thus  the  WT  would 
be  expected  to  be  equally  effective  for  various  transmis¬ 
sion  geometry  patterns  and  frequencies,  if  the  sound  speed 
profile  is  constant  and  the  bottom  fully  absorptive.  There 
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SOURCE  OEPTH  DETERMINATION  CURVE 
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Figure  26.  vf Detsraiaation  Curve;  F 
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SOURCE  OEPTH  OETERM I NATION  CURVE 
ISOSPEED  FREQUENCY  VARIATION 


Figure  27.  Iso^peed  Frofile  SD  Determination  Curve 
Variation  -2. 
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Figure  28.  Iso^pefd  Profile  SD  Deteroinaticn  Curve:  HD 


■ay  be  slight  effects  for  receiver  depth  and  frequency,  but 
these  have  not  been  substantiated. 

C.  BEALISTIC  PBOFILE  CASES 

1.  SEaiiAl  laiialiEs 

Fig.  29  is  the  conpatison  of  spatial  variability  for 
the  curves  which  resulted  from  the  spectra  of  profiles  A  and 
B,  (The  sound  speed  profiles  correspoading  to  the  pre¬ 
diction  curves  will  be  included  for  each  comparison.)  These 
were  both  sampled  on  November  16,  1980,  and  represent  a  spa¬ 
tial  separation  of  about  200  nm.  Both  profiles  result  in  a 
slight  bias  below  the  theoretical  Lloyd’s  Mirror  curve. 
There  is  more  variation  along  each  curve,  which  is  to  be 
expected  because  the  profiles  are  not  isospeed.  Also,  par¬ 
tially  absorbing  bottoms  are  now  used  as  compared  to  the 
fully  absorbing  bottoms  used  '^ith  the  isospeed  cases.  The 
profile  and  bottom  effects  have  not  yet  been  separated  to 
determine  which  contributes  more  to  the  departure  from 
theory.  This  would  require  an  extensive  set  of  intercompari¬ 
sons  between  these  two  parameters.  (The  present  intent  is 
to  show  that  ocean  variability  is  important.)  It  is  diffi¬ 
cult  to  differentiate  between  the  two  curves,  but  it  appears 
that  there  is  a  slight  trend  in  B  toward  smaller  null 
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spacings.  Sound  speed  profile  B  is  slower  than  K  over  the 
entire  depth.  A  slower  sound  speed  would  appear  to  indicate 
that  a  ray  of  lesser  angle  would  be  required  to  provide  the 
second  arrival  for  the  interference  pattern.  That  would  cor¬ 
respond  to  a  greater  vertical  wavenuaber  and  decreased  beta 
null  spacing. 

Fig.  30  compares  profiles  C  and  D,  which  were  sam¬ 
pled  at  the  sane  locations  as  profiles  A  and  B,  respec¬ 
tively,  but  on  December  6,  1980.  A  similar  pattern  is  noted 
with  D  spacing  less  than  C,  which  may  again  correspond  to 
slower  sound  speeds  fcr  profile  D. 

2.  Temporal  Vasiatjon 

Fig.  31  compares  profiles  A  and  C  which  represent 
the  same  location,  but  a  separation  of  20  days.  There 
appears  to  be  a  curve  crossing  at  about  800  feet,  such  that 
spacing  is  less  at  A  for  sha? .  depths  and  greater  for  A 
below  800  feet.  The  A  sound  spped  profile  is  slower  than  the 
C  profile  down  to  800  feet  and  the  pattern  is  reversed  from 
800  to  2000  feet.  Slower  speeds  appear  to  be  related  to 
shorter  null  spacings,  which  appears  consistent  with  this 
example.  The  variaticns  between  the  curves,  however,  are 
quite  small. 
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Figure  30.  Spatial  Variation  (C  ,0) :  SO  Dsteroination  and 
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Fig.  32  presents  a  sieilar  coeparison  for  profiles  B 
and  D.  The  speed  profiles  show  aleost  no  differences,  except 
that  D  is  slightly  slower.  This  result  appears  as  a 
slightly  sharper  deep  sound  channel  which  will  mean  slightly 
slower  welocity  for  the  direct  path.  That  effect  may  be 
reflected  as  a  slightly  greater  null  spacing  as  shown. 

3.  BgL^ge 

Results  for  profile  A  war*^  compared  for  ranges  of  50 
and  100  nm  in  Fig.  33.  Mo  significant  variations  were 

noted,  which  is  similar  to  the  results  for  *h6  isospeed 
cases. 

4.  Receiver  Depth  Variation 

Results  for  profile  A  were  also  compared  for 

receiver  depths  of  300,  1000  and  9000  feet  in  Fig.  34. 

There  may  be  slightly  better  agreement  for  the  shallow 
receiver,  but  this  is  not  certain. 

D.  COMPOSITE  COMPARISON 

Fig.  35  is  a  composite  of  all  isospeed  and  actual 

profile  data  sets.  Isospeed  data  show  slightly  better 
agreement  with  theory;  however  realistic  profiles  are  a  rea¬ 
sonable  approximation  to  the  theoretical  curve.  Realistic 
profile  variations  from  the  isospeed  curve  appear  to  be 
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related  to  the  sound  speed  profile  and  different  bottom 
boundary  conditions.  The  general  relationships  of  the  data 
indicate  that  below  about  800  feet  no  significant  variations 
are  noted.  Specific  curves  will  apparently  provide  better 
predictions  at  shallower  depths.  It  may  be  sufficient  for 
certain  applications  to  provide  climatological  or  theoreti¬ 
cal  predictlcn  curves.  These  curves  con  Id  be  enhanced  or 
modified  by  updated  curves  that  would  reflect  the  current 
sound  speed  profile. 
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Figure  33.  Range  (R)  Variation  for  Profile  A 
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Figure  35.  Composite  Comparison  for  Isospsed  and  Realistic 
Cases. 
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V.  COWCLg Slows 


The  WavenaBber  Technique  (HT)»  as  proposed  by  Lauer 
(1979)  has  been  investigated  for  source  depth  determination. 
Several  idealised  and  actual  scenarios  of  oceanic  and  acous¬ 
tic  parametric  variations  vere  tested.  Ihere  is  excellent 
agreement  between  the  ideal  test  cases  and  the  theoretical 
Lloyd's  Mirror  effect.  The  realistic  cases  showed  signifi¬ 
cant  similarity  as  well,  but  with  some  bias  toward  smaller 
wavenumber  spacing.  This  is  gualitatively  attributed  to 
sound  speed  profile  variations  in  the  water  and  sediment. 

There  appears  to  be  a  response  of  the  WT  to  variations 
in  input,  which  is  beneficial  in  allowing  evaluation  of  the 
relative  importance  of  individual  parameters.  Changes  in 
range  and  water  column  depth  do  not  indicate  any  significant 
variation.  Some  differences  are  possibly  suggested  for 

receiver  depth  and  frequency,  but  the  results  have  not  been 

• 

fully  analyzed.  The  WT  appears  to  respond  well  to  changes 
in  the  environmental  parameters,  such  as  sound  speed 
profile.  This  is  seen  from  the  average  bias  of  the  actual 
cases  from  the  isospeed  cases.  The  bias  does  appear  to  be 
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relatlvaly  consistent  and  shoes  the  importance  of  oceanic 
▼ariations  in  both  space  and  in  tine.  Thus,  it  appears  that, 
if  information  is  available  for  the  bottom  structure,  and  if 
a  realistic  sound  speed  profile  can  be  generated,  then  a 
very  specific  family  cf  curves  could  be  generated  for  vari¬ 
ous  frequency  and  receiver  depth  combinations.  This  simplic¬ 
ity  indicates  possible  operational  applicability  and 
flexibility. 

These  are  the  first  known  source  depth  determination 
curves  for  the  WT  and  thus  there  is  no  comparison  available. 
It  nay  be  that  there  are  some  biases  in  these  curves  which 
would  be  discovered  when  comparative  runs  are  made  using 
other  inputs  or  propagation  models.  Comparison  with  observed 
signals  might  also  reveal  biases  that  are  related  to  the  'AT 
or  to  the  propagation  model.  An  explanation  has  not  been 
determined  for  the  (I- shape  nature  of  the  wavenumber  spectrum 
plots  and  thus  slightly  different  curves  might  result  if 
this  is  found  to  bias  the  null  spacing.  Also  there  is  a 
bias  associated  with  the  method  of  manually  estimating  the 
null  spacing.  This  might  be  resolved  with  a  computer  proce¬ 
dure  to  locate  the  nulls.  The  experimental  curves  actually 
agree  well  with  theory  under  the  assumption  of  a  flat. 
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totally-reflecting  surface  and  a  fully-absorbing  bottoa. 
The  premise  that  the  theory  night  be  approximately  accurate 
in  realistic  cases  appears  valid.  This  overall  ganerally 
snooth  fit  supports  the  validity  of  the  wr. 

This  has  been  an  initial  investigation  of  the  HT,  and 
thus  only  a  broadbrush  study  was  possible.  These  prelimi¬ 
nary  results  indicate  that  additional  analysis  might  be  ben¬ 
eficial.  The  response  of  the  BT  to  input  from  different 
propagation  models  would  determine  if  results  from  different 
models  could  be  compared.  Detailed  investigation  of  the  spe¬ 
cific  effects  of  various  parameters  would  lead  to  a  quanti¬ 
tative  analysis  of  the  sensitivity  of  the  WT.  An  explanation 
for  the  anomalous  O-^aped  spectrum  and  the  development  of 
an  algorithm  which  can  accurately  determine  null  locations 
would  be  necessary  before  the  HT  could  be  used  practically. 
An  extremely  important  investigation  would  be  the  comparison 
of  analyzed  actual  signals  to  predictions  from  the  MT  to 
determine  applicability,  accuracy,  consistency  and  any 
appropriate  biases. 

The  intent  of  this  research  was  to  determine  the  valid¬ 
ity  of  the  BT  by  development  of  source  depth  deterninaricn 
curves  and  to  evaluate  the  oceanic  and  acoustic 


sensitivities  of  these  curves.  The  results  obtained  indicate 
that  in  general  the  UT  appears  to  be  valid  wi^h  respect  to 
the  Lloyd’s  Mirror  effect.  The  parasetric  variations  which 
were  investigated  show  that  oceanic  variability  is  ertreeely 
isportant,  but  that  the  «T  is  flexible  for  other  paranetars, 
such  as  transiission  geoeetry.  The  »T  appears  to  be  viable 
for  operational  applications.  Phe  additional  investigations 
suggested  above  would  deteraina  the  validity  of  this  suppo¬ 
sition.  It  is  conceivable  that  the  WT  could  be  developed 
into  an  operational  prediction  product. 
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